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The weight distribution of the code vectors of a binary group code has 
been referred to as the spectrum of the code. This paper presents a technique 
for calculating the spectrum of such a code, the spectra of shortened codes 
obtainable from the code, and what are defined as the level weight structures 
of the code. 

The method is conceptually straightforward and readily adaptable to 
digital computers. It involves operations no more complex than the addition 
of two (n — k) -tuples, the determination of the weight of certain (n — k)- 
tuples, and the ordinary addition of certain integers. Its computational 
complexities are independent of the code parameters. In principle, it may 
be used for any binary group code, but it is particularly useful for codes in 
which the number of parity check positions per code vector is rather small 
although the number of information positions may be large. 

I. INTRODUCTION 

The need for reliable data transmission systems has prompted the 
investigation of various coding techniques which attempt to detect 
and/or correct transmission errors. Because of the relative ease with 
which binary codes can be implemented, these codes have received 
special attention. It is with certain properties of these codes that this 
paper is concerned. 

In general, the encoder receives a block of k binary symbols (called 
a message) from a message source from which it determines (n — k) 
binary parity check symbols (called an ending). The message symbols 
and the ending symbols may be interleaved or transmitted sequentially 
thus forming a block of length n (called a code vector). Because any 

* The material presented in this paper formed Appendix II of the dissertation 
"Coding for Numerical Data Transmission" submitted by the author to The Johns 
Hopkins University in conformity with the requirements for the degree Doctor of 
Philosophy. 
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code in which these symbols are interleaved is equivalent 1 to a code in 
which the message and ending are transmitted sequentially, attention 
may be restricted to the latter situation. 

The elements and 1 form a field. Two vectors (or n-tuples) whose 
components are these field elements may be added by adding modulo 
2 the corresponding components of each vector. The symbol © will 
be used to denote this addition of vectors. 

The set of all possible n-tuples forms a vector space V» of dimension 
n over the field of two elements. A subset V is said to form a group 
code if the n-tuples in the subset form a group. Over the field of two 
elements, a set of vectors that forms a group is a subspace of V n . There- 
fore, the vectors of a group code form a subspace of V n . 

The weight of a vector u is the number of nonzero components in u 
and is denoted by w[u). The distance 1 between two code vectors u and 
v is w[u © v]. Because the code vectors form a group, there exists a code 
vector t = u © v. The distance between u and v is thus equal to w[t]. 

Because of this relationship between code vector weights and distances 
between code vectors, it is useful in evaluating the error detecting and/or 
correcting capabilities of group codes to be able to determine the num- 
ber of code vectors of each possible weight — i.e., from to n. This 
information has been called the spectrum of a code and can in principle 
be obtained by calculating in detail each of the possible 2 code vectors 
and then determining the weight of each of these code vectors. However, 
this method is not computationally feasible for values of k which are 
most often of interest. 

MacWilliams 2 has determined a system of linear equations which 
relate the set of integers that forms the spectrum of a given code to the 
set of integers that forms the spectrum of its dual code. The method is 
particularly effective for codes in which the dimension of the dual code 
is relatively small so that the spectrum of the dual code is readily ob- 
tained. 

The method presented herein enables the direct computation of the 
spectrum of a code without the actual formation of every code vector. 
The technique also gives both the spectrum of each of the possible 
shortened codes which may be obtained from the given code and the 
level weight structures of the given code. The level weight structures 
(which are defined later in this paper) have proved useful in the study 
of the effectiveness of error-correcting codes for numerical data trans- 
mission 3 and may indeed be of interest in other areas of code evaluation. 

The method is conceptually quite simple, readily implemented on a 
digital computer, and does not depend upon the solution of any equa- 
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tions. In fact, the only operations involved are the component by com- 
ponent modulo 2 addition of (n — A:) -tuples, the determination of the 
weight of certain (n — fc)-tuples and the ordinary addition of certain 
integers. 

II. COMPUTATIONAL TECHNIQUE 

Let k denote the dimension of the code space V and let Ej (1 ^ j ^ k) 
denote the k basis vectors of V. Take Ej in the usual systematic form 

Ej = e j \C j (1) 

where the message e,- is the fc-tuple with a 1 in position j and all other 
positions and Cj is the (n — A:) -tuple ending assigned to the message 
ej . Note that if the code is specified by a parity check matrix 1 in the 
form 

H = (hh • ■ ■ hl n . k ) (2) 

where hi (1 ^ i ^ k) is the column of H in the ith position and /„_* 
is the (n — k) X (n — k) identity matrix, then Cj is simply the trans- 
pose of hj and Ej is readily obtainable. 

The vectors E\ , E% , • • • , E\ generate a subspace of V of dimension I 
which we shall denote as T t . Y k is the code itself and T t is the set of 
code vectors in which information positions / + 1, / + 2, • • • , A" are 0. 
r is defined as consisting exclusively of the all code vector. 

Let A; = Ti — Ti-i-Ai , which is called the /-level of the code, is the 
set of code vectors in which information positions I -4- 1, 1 + 2, • • • , fc 
are and information position 2 is 1. Any code vector in A* is the sum of 
Ei and some vector in I\_i . 

The basic idea is to form for Ti an ending -weight matrix S . For 
convenience we shall deviate from usual practice and number the rows 
and columns of S ( ' beginning with 0. The entry s a ,^ (0 in row a and 
column /3 of S {1) denotes the number of code vectors in Ti of weight a 
whose endings are the (n — A) -bit binary representation of /3 (denoted 
by B(/3)). There must be (n + 1) rows in S U) to allow for all possible 
code vector weights and 2"~ columns to allow for all of the possible 
(n - A;) -bit endings. Therefore, S U) is an (n + 1) X 2"~ k matrix. 

The utility of this technique lies in the ease with which S (l) may be 
obtained from S { ~ u . Suppose that $ l -1 is known. The code vectors of 
A; are formed by adding Ei to the code vectors of T;_i . However, the 
special form of Ei makes this operation equivalent to placing a 1 in 
information position I of each vector in T/_i and, at the same time, 
adding C t to the ending of each code vector in I\_i . 
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Any code vector of weight a in IVi whose ending is B((3) becomes a 
vector in A/ with ending B (/3) © Ci and of weight 7 where 

7 = a + 1 + w[B(p) © Ci] - w[B(j3)]. (3) 

For those values of a (0 ^ a ^ n) and 7 (0 ^7 ^ n) for which (3) 
may be satisfied, 



S 7 ,B-1[B(/S) e c,] 



— S 7 ,B-i[B(/j) e C|] T- S a ,/s 



where ZT 1 (the inverse of B) is the operator such that \p = B~ B\$\. 

In general, it is not possible to satisfy (3 ) for every 7 (0 ^ 7 ^ n ) . 
However, because all code vectors in IVi whose endings are B(fi) (i.e., 
the code vectors giving rise to the nonzero entries in column /3 of S -1 ) 
become code vectors in A* of weight in the range through n, all values 
of a corresponding to nonzero entries in column /3 of S produce 
values of 7 such that ^ 7 ^ n. For these values of 7, (4) may be 
applied. 

On the other hand, values of 7 which would require values of a outside 
of the range ^ a ^ n in order to satisfy (3 ) are those values of 7 for 
which it is impossible to have code vectors in IVi of weight a whose end- 
ings are B(f3). For these values of 7, 

(I) _ d-l) /c\ 

S 7 ,fl-l[B(0) e C,] — S 7 ,B-1[B(0)©C/] • \<J) 

The column numbers referred to in (4) and (5) are independent of 
a. Furthermore, as a increases, (4) and (5) simply refer to different 
elements in the same column. For this reason, these results may be 
expressed as column operations thus leading to a conceptually simple 
result. 

Let Sfl ( ° denote column /3 in S {1) . Define <r 0) to be a shifting operator 
which, when applied to sp l) , shifts each element of sp 1 ' by j positions 
filling in any resulting blank positions with zeros. For example, if 

"0" 






then 



(2) „ U) _ 

a -Sp = 
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and 






In terms of a , the relationships expressed in (3), (4), and (5) may be 
conveniently expressed as 

(J) _ „ C-l) i («>[B(0)ffiC,]-i0[B<fl)]+l) „ (l-l) /«v 

Clearly all columns of S (l) are obtained by successively applying (6) 
as |9 runs from to 2 n_ * — 1. 

It is important to notice the great simplicity of (6). In practice, it 
involves shifting one column of S _1 and then combining by simple 
addition the elements of this column with those of another column of 
<S ( ' _1) to obtain a column of S \ Determining the number of positions 
that as should be shifted and the column with which s/j ( ' -1) should 
be combined is extremely easy. In particular, the operations in (6) are 
readily adapted to digital computer operations. 

If S (0) is known, the remaining ending-weight matrices can be succes- 
sively obtained. The only code vector in To is the all vector. 
Therefore, 

e (0) — i 
and 

8.J* = 

for all other values of a and 0. 

Now that the method for constructing the ending-weight matrices 
has been presented, the following will serve to indicate how the desired 
information is extracted. 

(?) Spectrum of the code: The total number of code vectors of 
weight a in the code is 

"if «M W . (7) 

0=0 

The spectrum of the code is obtained from S lk) by using (7) 
for each value of a(0 ^ a ^ n). 
(ii) Spectra of shortened codes: Let k' denote the number of in- 
formation positions in the shortened code — i.e., k — k' in- 
formation positions are deleted. Assume that the deleted posi- 
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tions are information positions fc + 1, fc + 2, • • • , fc. The 

total number of code vectors of weight a in the shortened code is 

2n-i_i 

S ij". (8) 

(3=0 

The spectrum of the shortened code is obtained from S lk,) by 
using (8) for each value of a(0 ^ a ^ n — k + fc ). 
(in) Level weight structure: The set of code vectors A* has been 
referred to as the Mevel code vectors and the weight distribution 
of these code vectors as the Z-level weight structure. Note that 
the Mevel weight structure is the difference between the spec- 
trum of the shortened code consisting of I information positions 
and the spectrum of the shortened code consisting of (I — 1) 
information positions. 

Let ni, a denote the number of code vectors of weight a on 
the Mevel. The number of code vectors of weight a in Ti-i is 

2n-k— i 

E (1-1) 

S a>/ 3 
(3=0 

Similarly, the number of code vectors of weight a in Ti is 

2n-fc— i 

s a,0 • 
(3=0 

It follows that 

2 n-*_i 2 n ~ k — 1 

n l>a = E s a J l) - £ a*/" 1 '. (9) 

0=0 (3=0 

III. CONCLUSIONS 

The spectrum of any group code, the spectrum of any shortened code, 
and all level weight structures are obtainable in a straightforward 
manner by means of operations no more complex than the addition of 
two (n — k) -tuples (to determine the columns to combine), the com- 
putation of the weight of certain (w — A) -tuples, and the repeated 
addition of integers two at a time (to actually combine the columns). 
The number of computations does depend upon the parameters n and k 
but the method has the advantage that the complexity of the operations 
is invariant. Because the number of computations and the number of 
computer storage locations required for the ending-weight matrix are 
sensitive to changes in (n — k) but rather insensitive to changes in fc, 
the method is most effective for codes in which in — k) is moderate al- 
though k may be quite large. 
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As presented, the method treats each of the k ending-weight matrices 
in a similar manner by determining all of the (n + l)-2" _fc entries of 
each matrix. Computing time can be saved by realizing that the maxi- 
mum possible weight of a vector in Ti is I + n — k and, thus, that it is 
only necessary to compute the first I + n — k rows of S because the 
remaining rows contain zero entries exclusively. Additional programming 
sophistications, including processing only those columns of *S ( _1) which 
contain nonzero entries in obtaining *S ( ' (particularly for the smaller 
values of I), improve the computing efficiency of the method. 

This technique was originally developed for computing the level 
weight structures of certain codes. Thus, if the level weight structures 
and/or the spectra of the shortened codes are desired, this method offers 
a straightforward and effective means of obtaining such information. 
However, if all that is desired is the spectrum of the code, then under 
some conditions the method developed by Mac Williams 2 may be prefer- 
able from a computing time point of view although the conceptual sim- 
plicity of this method is still appealing. In any case, the relative ad- 
vantages of the two methods should be considered before deciding which 
to use for a specific application. 

The method has been used successfully to compute the level weight 
structures and the spectra of the (15,11), (31,26), and (63,57) Hamming 
perfect single error-correcting codes. In each case the information was 
obtained on an IBM 7094 digital computer in less than 0.01 hours. 

IV. NUMERICAL EXAMPLE 

The parity check matrix for a (7,4) Hamming code is 

/l 1 1 1 0\ 
H = 1 1 1 1 . 
\1 1 1 1/ 

The basis vectors for this code are 

Ei = 1000 111 

E 2 = 0100 110 

E 3 = 0010 101 

Ei = 0001 011 

message ending. 

The sets of code vectors referred to as T 3 and A 4 are listed in Table I . 
The appropriate values of a and are given next to each vector. 
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Table I 






Ti 


a 





Al 


a 





0000 000 








oooi on 


3 


3 


1000 111 


4 


7 


1001 100 


3 


4 


0100 110 


3 


6 


0101 101 


4 


5 


1100 001 


3 


1 


1101 010 


4 


2 


0010 101 


3 


5 


0011 110 


4 


6 


1010 010 


3 


2 


1011 001 


4 


1 


oiio on 


4 


3 


0111 000 


3 





1110 100 


4 


4 


1111 111 


7 


7 



Tabulating this information yields *S (3) and S w . 



S< 3 >:« 



12 3 4 5 6 7 



12 3 4 5 6 7 



10 
00000000 
00000000 
110 110 
00011001 
00000000 
00000000 
00000000 



S™:< 



10 
00000000 
00000000 
11111110 
1111111 
00000000 
00000000 
1 



We now turn to use the method herein developed to obtain S w 
from S {3) . Specifically, we use (6) first with = and then successively 
increase /3 until fi = 7. 

When = 0, £ _, [£(0) © C 4 ] = 3. Thus, (6) reduces to 



S 3 



(4) 



— o (3) 4- ,r (3) c (3) 
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Lo_ 




uO 



which is indeed correct. 

Now let/8 = 1. Then E~\B{\) C 4 ] = 2 so (6) yields 






(3) 



■ (1) „ (3) 

+ a -Si 
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+ 
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l(L 
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The remaining columns of S w are obtained in a like manner as /3 in- 
creases to 7. 

The spectrum of the code is obtained by summing across the rows 
of S w . The spectrum of the shortened code resulting from the deletion 
of the fourth information position is obtained by summing across the 
rows of *S (3) . The 4-level weight structure is the difference between these 
spectra. This information is tabulated in Table II. 





Table II 




Weight 


Code Spectrum 


Shortened Code 
Spectrum 


4-Level Weight 
Structure 





1 


1 





1 











2 
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7 


4 


3 


4 


7 


3 


4 


5 











6 











7 


1 





1 



For an illustrative example, it was necessary to confine ourselves to 
a code in which k is small. However, it should be realized that the true 
utility of the method lies in the fact that it can, without modification 
or additional complexity, be used for codes in which k is quite large. 
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